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We demonstrate that Bragg scattering from a condensate can be sensitive to the spatial phase 
distribution of the initial state. This allows preferential scattering from a selected spatial region, 
and provides a robust signature for a vortex state. We develop an analytic model which accurately 
describes this phenomenon and we give quantitative predictions for current experimental conditions. 
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Bragg scattering of matter waves from laser beams 
is a well known phenomenon in atom optics [e.g. 
and has proved to be an effective tool for manipulat- 
ing and analysing Bose condensates. In that context it 
has been used in experiments to coherently split a con- 
densate 1^,^, and to perform momentum and phonon 
spectroscopy Theoretical calculations show it may 

be used to extract quasiparticle amplitudes from weakly 
excited condensates Two groups have constructed 
Mach-Zehnder type interferometers using sequences 
of Bragg pulses to probe condensate phase |p|-pT[| . In all 
of those cases the Bragg scattering process is most eas- 
ily understood in terms of a momentum space process in 
which the Bragg fields transfer momentum ftq to each 
atom. In a beam-splitter, half the initial amplitude is 
translated in momentum space by ftq, while in momen- 
tum spectroscopy a selected narrow group of momentum 
components is translated by hc{. The spatial structure 
of the initial condensate plays no particular role in the 
Bragg process for those cases. 

A major interest in condensates however, is that their 
macroscopic wavefunctions may be engineered into spe- 
cific spatial structures. For example, the vortex state, 
which has been recently observed by two groups |]l2| , |l3| , 
has a phase circulation of 27r about a vortex core. The 
most direct means for observing such a phase distribu- 
tion is by interference with a separate well characterised 
matter field [e.g. |14 1^], but this may not always be 
experimentally convenient. In this paper we show that 
under appropriate conditions Bragg scattering is sensi- 
tive to the spatial phase distribution of the initial con- 
densate and therefore allows preferential scattering from 
a selected spatial region. In the case of a vortex, this 
gives a distinctive signature which we illustrate in Fig. 
|l|. There the density distribution of a trapped conden- 
sate vortex state is shown following application of Bragg 
pulse chosen according to criteria developed later in this 
paper. In Fig. [^(a) the streaming output emerges from 
only one side of the initial condensate, and gives an asym- 
metric density pattern which should be detectable with 
current experimental technology. In Fig. [l](b) where the 
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FIG. 1. Spatial density of a 2D vortex state at time 
t = Q.&tf) after excitation by a Bragg pulse with frequency 
detuning (a) uj = 21&ujt and (b) uj = 203tJT (i^t is the har- 
monic trapping frequency). Other parameters are V = G.2u}t, 
q = 14/2:0 and w = SOOuiq. Units: time to = 1/<^t; distance 
xq = h/2mujT; coUisional interaction strength wo = hujTXQ. 
The dashed line denotes a suitable position for measuring the 



output beam profile (see text). 



frequency difference of the laser fields has been changed, 
the streaming output field has a density node that is an 
order of magnitude wider than the vortex core itself (i.e. 
the healing length in the centre of the original conden- 
sate). In the following we develop a treatment of such 
spatially selective Bragg scattering which explains the be- 
haviour in terms of the underlying spatial phase sensitive 
mechanism and we give analytic solutions appropriate to 
3D condensates. We treat the process of Bragg scat- 
tering from a condensate using a previously developed 
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formalism , in which the condensate mean- field ^ (in 
an interaction picture) evolves in a far detuned light-field 
grating formed by two plane- wave laser beams, according 
to the equation 



+ VT{r) 

2m 

hV{t) cos(q • r - wt)^'. 



(1) 



Here Vt{v) is the trap potential (assumed to be harmonic 
in the examples of this paper) and w ( = ANirh^a/m) is 
the usual nonlinearity parameter of the Gross-Pitaevskii 
equation. The final term in Eq. (|l]) describes the inter- 
action with the Bragg grating: lo and ftq are the two- 
photon detuning and recoil momentum respectively and 
V{t) = |0(i)p/2A is twice the ac-stark shift at the inten- 
sity peaks. For the scattering to be in the Bragg regime 
the pulse length Tp of the lasers must be long enough 
that the Bragg resonance is resolved (i.e. Tp 3> : 
where ojq = fuf' /2m is the recoil frequency). In order to 
obtain spatially selective Bragg scattering it is also neces- 
sary that the average value of V is sufficiently small that 
the total amount of condensate scattered in time Tp is 
small compared to the unscattered amount. In addition 
we assume the Bragg pulse length (and the time of ob- 
servation) is shorter than a quarter trap period, to avoid 
the trap forces significantly altering the momentum of 
the scattered beam. 

Equation (Q) can be directly numerically simulated, as 
for example we have shown in the two dimensional case of 
Fig. |l|. We have also obtained an analytic solution, which 
provides detailed insight into the underlying physics and 
enables us to provide quantitative calculations over a 
wide range of parameters in the three dimensional case. 
For the analytic solution, we assume that the recoil mo- 
mentum h(\ of the Bragg grating is much larger than the 
momentum width of the initial condensate (which is cen- 
tred about zero momentum), and also assume that the 
detuning u is close to the Bragg resonance. This means 
the scattered wavepacket is well separated in momentum 
space from the initial state and a slowly varying envelope 
approximation can be made mM so that we write 



vE'(r,t) = V'o(r,t)+^Ai(r,t) e^^'i--"*). 



(2) 



In Eq. (H) Tpi is the scattered wavepacket while i/'Oi which 
we call the mother condensate, represents a condensate 
with a momentum wavepacket centered on zero. At i = 0, 
"00 is exactly the initial condensate. Substituting Eq. 
(^ into Eq. ([^) and projecting into orthogonal regions 
of momentum space, we obtain to first order in ipi the 
coupled equations 
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where Q = hq/m is the velocity of the scattered atoms, 
S = uj — ujq and we have neglected the rapidly oscillat- 
ing term 'ipo4'i- Our interest is in the scattered state 
■01 which is assumed to be small, so we can neglect the 
scattering from 0i back to 0o (i.e. V{t)^i in Eq. (||)). 
This means that 0o then evolves according to the usual 
Gross-Pitaevskii equation, and typically we choose 0o to 
be an eigenstate of that equation. We can also ignore the 
term in Eq. (^ which describes momentum diffusion 
about the centre of momentum q of the ^'i packet, and 
is small on the time scales we consider. Eq. 
be solved to give 



can now 



0i(r,t) = I dse-'^^('-'*'^V(s)^o (r+Q(s-0,s) , 



(5) 



where 



K{r,t,s) ^ {s - t)S + - I ds' Vrir + Q{s' - t)) 
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Eq. (|5|) allows us to visualise the formation of the scat- 
tered state as follows. As the scattered packet moves 
across the mother condensate, the amplitude 0i at a 
given point (stationary in the frame moving with veloc- 
ity Q ) is built up from the sum of contributions coupled 
in from successive points along the mother condensate. 
The contribution coupled in at time s from a position 
r + Q(s — t) on the mother condensate evolves to time 
t with the propagator exp{—iK{r,t,s)). The scattered 
state at r and t will be appreciable only if the contribu- 
tions constructively interfere. Writing 0o(r,s) in terms 
of amplitude and phase Ao{r,s) exp{iSo{r, s)), then the 
condition for an appreciable scattered state to form is 
that the phase Q{r,t,s) = -K{r,t, s) + So{r + Cl{s' - 1)) 
be stationary. We take ipo to be a stationary eigenstate, 
i.e 0o(r, s) = Ao(r) exp(i(5o(r) — /ips)), where is the 
eigenvalue of the mother condensate. Furthermore we 
make use of the Thomas Fermi solution w[Ao]^ = fio + Vr 
so then 

K^{t-s)ifio-5) + j J' ds' [Ao{r + Q{s' -t) f . (7) 

The condition for stationary phase dQ{r,t,s)/ds ~ 
then gives the generalised Bragg resonance condition 



^[Ao(R)]' + Vr5o(R)-Q 



R=r+Q(s-t) 



(8) 



If O has sufficiently large spatial curvature, one particu- 
lar time s may dominate the stationary phase contribu- 
tion. However for the cases of a ground or vortex initial 
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state, & varies sufficiently slowly that most contributions 
from along the line R = r + Q(s — i) on the initial con- 
densate can be in phase. For an initial ground state, 
Eq. reduces in the linear case to the usual Bragg 
resonance condition 6 = 0, while for the nonlinear case, 
S = {w|V'o(R')P/^} ill which the braces indicate that the 
precise value of the shift is obtained by a suitable aver- 
age along the path R [p7[ . For a general initial state the 
second term in Eq. (fel), which can be interpreted as a 
Dopplcr condition, must be included. 

The full spatial solution for the scattered state con- 
tains a great deal of information, but we can extract a 
useful and compact signature by considering the steady- 
state density profile of the scattered beam at the edge of 
the initial condensate. In the 2D simulation of Fig. |l| for 
example, the steady-state density profile could be mea- 
sured along the dashed line. For definiteness, we shall 
take the case where q is in the y direction. We define the 
steady-state density profile of the scattered state to be 

Vq{x,z,uj) = \ipi{x,y = R,z, Ts)\'^, (9) 

where the distance R is sufficiently large that the density 
of the mother condensate in the plane {x,y = R,z) is 
negligible (i.e Vg « l^'P). For strongly interacting con- 
densates it is suitable to use R = Rtf, where Rtf is 
the Thomas- Fermi radius (in the xy-plane). rg, the time 
to reach steady state can be estimated as 2R/Q, i.e the 
time for the scattered state to traverse the mother con- 
densate. Since experiments would usually measure the 
density projected along the z axis (the line of viewing), 
we define the projected steady-state density profile 

Dq{x,Lu) = J dz'Dq{x, Z,Uj), (10) 

to characterise the results of the Bragg scattering. 

We begin by considering the case of two spatial dimen- 
sions, which contains the essence of the physics. Fig. || 
shows Dq{x, oj) calculated from the full 2D numerical so- 
lution of Eq. (|l|) , for cases where the initial state is (a) 
a noninteracting ground state (b) a condensate ground 
state, (c) a condensate vortex state. In Fig. || (d)-(f) 
we provide the comparison to the analytic solutions for 
the same cases, and it is apparent that agreement is very 
good. We note that in the analytic solutions a Thomas 
Fermi approximation is made for t/jq. For vortex states, 
the additional centrifugal potential m^/(a;^ + U^) is in- 
cluded and the phase is S'o(r) = m arctan(?;/a;). 

The results for the noninteracting ground state (Fig. 
^(a) or (d)) show, as expected, that the scattered beam is 
greatest (i.e. Bragg scattering is resonant) when oj = ujq 
(i.e. 6 = 0). The frequency width of Dq{x,uj) can be es- 
timated using a simple Fourier analysis on the integral in 
Eq. (|), to be Az5 - Q/Rtf- In Fig. |(b) [or (e)] the ef- 
fects of the condensate nonlinearity appear. At the centre 
of the scattered profile {x = 0) the resonant frequency for 
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FIG. 2. Steady-state Bragg scattered density profiles in 
2D for (a) noninteracting ground state (b) conden- 

sate ground state (w = SOOioo), (c) condensate vortex state 
{w = 500?i)o, rn = —1). Measurements are made at t = O.Gto 
for a Bragg field with V = IutjQ = 14/a;o. Frames (d)-(f) 
give the analytic solution of Eq. (Bl) for the cases correspond- 
ing to (a)-(c) respectively. Dotted lines, free particle resonant 
frequency {ujq = hq^ /2m); dashed lines, 2D nonlinear shifted 
frequency (a)„i = iOq + 2fio/3) @. 

Bragg scattering has been shifted by .84/i which is close 
to the value of {w|-0o(R)P/fi-} averaged along the cen- 
tre line of the mother condensate (=4/xo/5). The shift 
at the spatial edges of the scattered beam (x = ±Rtf) 
is less, because the mother condensate has lower average 
density along the appropriate lines R = r-t- Q(s — i), and 
so Dq(x,uj) has a crescent shape. The frequency width 
Ao (at a: = 0) is smaller than for the noninteracting 
ground state case (a) , due mainly to the increased spatial 
width Rtf of the mother condensate. The most signifi- 
cant result for this paper is Fig. ^(c), namely the vortex 
signature. At the 'resonant' frequency (indicated by the 
dashed line), the scattered density profile is essentially 
spatially symmetric (see also Fig. |^(b)), but at other 
frequencies the scattering is spatially asymmetric. This 
asymmetry, which we emphasize arises from the spatial 
phase asymmetry of the vortex, is robust, being present 
for a wide range of frequencies. The density node at the 
centre of the beam {x = 0) at the resonant frequency 
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also arises from the phase asymmetry: as the scattered 
wavepacket passes over the vortex core, the contributions 
from the mother condensate change phase sharply by tt 
and thus cancel. Our results can also be extended into 
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FIG. 3. Projected steady-state Bragg scattered density 
profiles from a condensate of 2x10^ Rb®'^ atoms in an m = — f 
central vortex in (a) prolate trap with uit — 2n x 50Hz, 
aspect ratio A = \/S, /iq ~ SI.Glut', (b) oblate trap with 
LOT = 2nxl00Hz, X = 1/10, /io ~ 9.5ujt- Dotted lines, free 
particle resonant frequency {ujq — 2ttx 15kHz); dashed lines, 
3D nonlinear shifted frequency {u>ni = + 4/xo/7). Bragg 
field {V — 2nx50Hz) provided by counter-propagating lasers 
of approximately 780nm. 

three dimensions by using the analytic solution Eq. (|^). 
In Fig. ^we present the behaviour of Dq{x, to) from a vor- 
tex condensate of 2x10^ Rubidium atoms in both oblate 
and prolate traps. The features discussed in the previ- 
ous paragraphs are unchanged in three dimensions. It is 
worth noting that our results can be readily connected to 
the well known dynamic structure factor S'(q, lo) which 
has been previously used to characterise the results of 
Bragg scattering [p|,p|Jl9|]. One can show, within pertur- 
bation theory, that 



dx Dq(x, Uj) . 



(11) 



In other words, the dynamic structure factor simply 
projects out all the spatial information of the steady- 
state density profile. 

Although we have concentrated on the steady state 
density profile our analytic solution contains other inter- 
esting results. For example, with a Bragg pulse of length 
T„ ^ Tc, the scattered state is itself a vortex, and a 



sequence of such pulses would produce a sequence of vor- 
tices streaming out from the mother condensate, until it 
becomes too depleted. 

In summary we have shown that under appropriate 
conditions Bragg scattering is sensitive to the spatial 
phase dependence of the initial matter field state, and 
therefore allows preferential scattering from a selected 
spatial region. We have developed an analytic model 
which accurately describes this phenomenon and explains 
the underlying mechanisms. When applied to a vortex 
state, a robust signature is obtained, for both oblate and 
prolate traps. 
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